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Abstract — This paper presents numerical simulations 
of dual-frequency second-order ultrasound field 
(SURF) reverberation suppression transmit-pulse 
complexes. Such propagation was previously studied 
in a homogeneous medium. Here instead the 
propagation path includes a strongly aberrating 
body-wall modeled by a sequence of delay-screens. 
The applied SURF transmit pulse complexes each 
consist of a high-frequency imaging 3.5 MHz pulse 
combined with a low-frequency 0.5 MHz sound speed 
manipulation pulse. Furthermore, the feasibility of 
two signal post-processing methods are investigated 
using the aberrated transmit SURF beams. These 
methods are previously shown to adjust the depth 
of maximum SURF reverberation suppression within 
a homogeneous medium. The request of the study 
arises because imaging situations where reverberation 
suppression is useful are also likely to produce pulse 
wave- front distortion (aberration). Such distortions 
could potentially produce time-delays that cancel the 
accumulated propagation time-delay needed for the 
SURF reverberation suppression technique. Results 
show that both the generation of synthetic SURF 
reverberation suppression imaging transmit-beams, 
and the following post-processing adjustments, are 
attainable even when a body-wall introduces time- 
delays which are larger than previously reported 
delays measured on human body-wall specimens. 

The peer-reviewed version of this paper is published 
in IEEE Transactions on Ultrasonics, Ferroelectrics 
and Frequency Control, vol. 59, no. 11, pp. 2588-2595, 
November 2012. DOI: 10 . 1109/TUFFC . 2012 . 2494 
The final version is available online at 
http://dx.doi.org/10.1109/TUFFC.2012.2494 The current 
e-print is typeset by the authors and differs in e.g. 
pagination and typographic detail. 



I. Introduction 

This work concerns dual-frequency transmit-beams that 
are utilized for reverberation noise suppression in ultra- 
sound image reconstruction by use of methods analyzed in 
[l], Q . The developments of the current paper regard the 
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feasibility of such beam generation within an inhomoge- 
neous medium that generates sound-speed variations that 
aberrate the beam. Although the subject is connected, 
here we don't apply or develop aberration correction 
algorithms. 

The quality of medical ultrasound images varies greatly 
between patients. Spatial inhomogeneities in compressibil- 
ity and density within the ultrasound propagation path 
cause additive reverberation (multiple scattering) noise, as 
well as resolution degradation due to wave-front distortion 
(aberration). This is especially prominent when imaging 
through ejr. the abdominal wall or the chest wall of obese 
patients [3|47j . A recent simulation study supports the hy- 
pothesis that the image degradation due to reverberation 
noise is more prevalent than degradation due to aberration 
noise HQ. 

In dual- frequency second-order ultrasound field (SURF) 
reverberation-suppression imaging, a synthetic transmit- 
beam is generated from the difference between two high- 
frequency (HF) imaging pulses transmitted in the same 
direction. This beam is in the following simply denoted 
the SURF beam. Both transmissions comprise a dual- 
frequency band pulse-complex with the HF part used 
for image reconstruction, added to a low-frequency (LF) 
part used for material compressibility manipulation. The 
high frequency is typically of the order 10 times the 
low frequency. The compressibility is pressure-dependent 
due to material nonlinearity. This makes the sound- 
speed pressure-dependent, being higher for the compressed 
medium than for the expanded medium flQ| . The LF 
polarity is switched for the second transmission, therefore 
making the LF pressure experienced by the second HF 
pulse opposite what was experienced by the first HF 
pulse. The pressure-dependent sound-speed thus causes 
the two HF pulses to propagate at different speeds. With 
increasing distance traveled, the difference in propagation- 
time required for the two HF pulses hence increases. 
Consequently, the difference obtained when subtracting 
the HF pulses is zero at zero depth and grows with 
covered distance to attain a maximum at 180° phase- 
shift. This is the effect that is exploited for reverberation 
suppression as follows. The SURF beam is the synthetic 
transmit-beam obtained when subtracting the propagated 
HF pulses at each spatial position. The amplitude of 
scattered (or reflected) pulses are severely reduced, hence 
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making the nonlinear sound-speed manipulation negligi- 
ble. Therefore scattered HF parts accumulate less relative 
time-shift than forward-propagating ones. Especially when 
the first scattering takes place at shallow depths, multiple 
scattered contributions within the received HF pulses thus 
give negligible contribution to the receive SURF difference 
HF field. The reverberation suppression ability of SURF 
imaging is illustrated by the fact that the beam is reduced 
near the transducer, as was studied in [1] where nonlin- 
ear numerical propagation simulations in a homogeneous 
medium were reported. In addition, the combination of 
tissue-harmonic pulse-inversion with SURF reverberation 
suppression was introduced. 

In case the dominating reverberation noise is due to a 
strong first scattering (reflection) taking place around a 
depth z a , the SURF synthetic transmit field may without 
modification or re-transmission of the sent-out pulses, 
be adjusted to increase the transmit-beam suppression 
around z a by use of post-processing methods as described 
and analyzed in [2| for a homogeneous medium. Both mul- 
tiple scattering and multiple reflection image artifacts are 
in the following considered as reverberation artifacts. The 
recent publications lH-[l3| concern SURF dual-frequency 
acoustics within a wider context. Dual frequency band 
pulse complexes have also been used for contrast agent 
detection (T34is(l - 

The purpose of the computer simulation study pre- 
sented in the current paper is to evaluate the feasibility 
of SURF transmit reverberation suppression beam gen- 
eration when a strongly aberrative body-wall is present 
within the wave-propagation path. Furthermore, it intends 
to compare these SURF beams to standard fundamental 
imaging transmit-beams, and in addition to test the post- 
processing SURF transmit field enhancement techniques 
described in [2[ in case of aberration. 

The study is important because the body-wall of a 
patient where standard fundamental ultrasound imaging 
is aggravated due to reverberation noise, is also likely to 
produce propagation time-shifts due to aberration. If these 
time-shifts heavily distort the accumulated time-shifts 
needed for adequate SURF synthetic beam generation, 
the reverberation suppression gain of the SURF imaging 
method is no longer obtained. 

This paper is organized as follows: First the dual- 
frequency SURF imaging and fundamental imaging exci- 
tation pulses used in the numerical simulations are pre- 
sented. Then a description of the body-wall model and 
field simulation method is given. The results section shows 
comparisons between transmit-beams from the described 
excitation pulses both for a homogeneous medium and 
after propagation through the body-wall. Beam profiles 
in the focal plane are derived from the transmit-beams, 
to quantify the effect the body-wall has on the transmit 
field. Finally, post-processing adjusted SURF transmit- 
beams are constructed for some chosen example depths 
of decreased transmit-beam amplitude. 



II. Theory and Methods 

A. Excitation pulses 

The two excitation pulses utilized to generate the SURF 
synthetic transmit-beam and the transducer apertures are 
equal as for a setup applied in [H, Q. The HF imaging 
part of each pulse complex has a center frequency at /h = 
3.5 MHz (50% fractional bandwidth at -6 dB) while the 
LF part, whose polarity is inverted for the second pulse, 
has a center frequency at /l = 0.5 MHz (25% bandwidth 
at — 6dB). The HF and LF surface pressure amplitudes 
are adjusted to keep the Mechanical Index (MI) below 1.9 
and are therefore set to 3.5 and 0.85 MPa respectively. The 
peak negative pressure value applied for MI calculation 
is read before filtering the SURF pulse complex, hence 
taking the amplitude contribution of the LF part into 
account. Therefore the highest MI is observed for the 
SURF complex of negative LF manipulation polarity. The 
axisymmetric aperture is focused at 82 mm with an outer 
radius for the HF part of 7.1mm and 10 mm for the LF 
part. The HF pulse transmission is delayed by To = 0.2 /is 
compared to the LF, measured at their centers. Fig. 2 in 
Ref. LU illustrates the influence of To on the transmit SURF 
pulse complex. 

The excitation pulse utilized to generate the comparison 
fundamental imaging field, is equal to the HF part of the 
SURF complexes described above. 

B. Body-wall model and pulse propagation simulation 

The applied forward wave-propagation simulation 
method solves the nonlinear Khokhlov-Zabolotskaya- 
Kuznetsov (KZK) wave equation taking attenuation and 
interaction between the HF and LF parts of the SURF 



pulse complexes into account [2CJ, |2l|. The transducer 
apertures modeled are rotational symmetric. Because a 
non-symmetric aberrating body- wall is taken into account 
within this work, all simulations are performed in three 
spatial dimensions, hence demanding longer computation 
times and greater storage capacity than for 2-D rotational 
symmetric simulations. 

The body-wall model utilized was generated using a 
set of 13 two-dimensional filtered time-delay white- noise 
screens, adjusted to emulate a strongly aberrating abdom- 
inal wall. The screens are equally spaced by Az = 3 mm 
thus giving a total wall thickness of d = 39 mm as 
shown in Fig. [TJ This body-wall is similarly modeled as 
in [22j], however the aberrative effect is here chosen to 
be more severe in order to test the feasibility of the 
SURF beam generation methods under harsher imaging 
conditions. The body-wall model construction method is 
thoroughly described in [23| and the characteristics of the 
modeled body-wall, as found from simulations of back- 
propagation of the signal from an excited point source 
in focus, are displayed in Table |L] At zero-pressure, the 
material property average values are: speed of sound 
1550 m/s, density 1.06 mg/mm 3 , and the nonlinearity 
parameter f3 n = 3.9. The attenuation is described by a 
power-law making it proportional to f 1A with the loss 
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Fig. 1. Body- wall model and simulation setup outline. The total wall 
thickness is d = 39 mm, the transmit focus is at F, and Az = 3 mm 
is the distance separating each delay screen. 

TABLE I. 

Characteristics of the simulated body-wall 



Amplitude correlation-length 1.4 mm 

Amplitude RMS fluctuation 4.5 dB 

Time-delay correlation-length 1.2 mm 

Time-delay RMS fluctuation 450 ns 



ao = 0.52dB/cm at 1MHz. Details on the simulation and 
attenuation implementations are given in Ref. [24]. 

Both the corresponding non-aberrated SURF beam and 
the non-aberrated fundamental beam are simulated and 
used as reference beams. All beams are generated by 
detection of the temporal maximum of the pressure field 
at each spatial coordinate. 

C. SURF signal post-processing 

Two different signal post-processing methods, further 
described and analyzed for a homogeneous medium in 
[H, are applied on the SURF field data filtered around 
the frequency /h = 3.5 MHz. They both aim to further 
suppress the SURF transmit-beam at some depth z a , so 
that the ratio between the beam energy within the imaging 
depth region and the energy at the chosen suppression 
depth is increased. This way image reverberation-artifacts 
due to first scattering from some strongly reflecting feature 
at the depth z a may be suppressed. 

In short, given that the two transmit HF waves prop- 
agated in conjunction with LF manipulation pulses of 
opposite polarities are s + (r,t) and s_(r, t), the second HF 
signal is processed to form the field s_(r,t) which is used 
to form the difference field 

SA(r,t) = s+(r,t) -s-(r,t). (1) 

The adjusted SURF synthetic transmit-beam is then gen- 
erated from sa(?j i). 

The first post-processing method involves a pure time- 
shift of the s-(f,t) field by r a so that 

§-(i f ,t)=s-(r,t-T a ). (2) 

The second post-processing method involves some more 
general filter \) Za which makes the on- axis pulse of s + and 
S- equal at the depth z a : 

s-(r,t) = t, za {s.(r,t)}. (3) 



Here this filter \) Za is determined by comparison of the 
two simulated fields as averaged over a laterally 0.5 mm 
wide region S on-axis at the depth z a and finding the 
filter that makes these two average signals equal. For the 
demonstrated inhomogeneous medium case, the adjust- 
ment filter \) Za utilized is not generated from the actual 
aberrated s + (r,t) and s_(r,t) fields, but instead from 
the s + (r,t) and s-(r,t) emerging from the homogeneous 
medium propagation. This approach is chosen to in some 
extent better emulate how an a priori determination of \) Za 
could be performed by experimental field characterization. 
The averaging over the region 5, which was not performed 
in [2[, is introduced to take a finite hydrophone width 
into account. The determination of \) Za in a real medical 
imaging setup may be more cumbersome, as further dis- 
cussed in [2[ . We remark that when receive channel-data is 
available, any of the proposed suppression methods may be 
applied, and images can be reconstructed independently 
for multiple choices of z a . Great flexibility in reverberation 
suppression customization is hence enabled. 

D. Compared transmit- beams and fields 

In the following is a listing of the simulated datasets 
utilized to generate the transmit-beams that are compared 
within this work: 

1) Non-adjusted SURF beam: The HF transmit fields 
s+(r,t) and s_(r, t), filtered around fn = 3.5 MHz on the 
same actual dataset as utilized in [l| , are subtracted with- 
out further modification to generate sa(^,£) = 5 +(r, — 
s-(r,i). 

2) Time- shift adjusted SURF beams: The time-shift 
adjustment post-processing method is applied on the same 
original dataset as used to generate the non- adjusted 
SURF, for a number of different z a to generate the ad- 
justed HF field s&(r,t) = s + (r,t) — s-(r,t — r a ). 

3) Filter- adjusted SURF beams: The filter- adjustment 
post-processing method is applied on same original dataset 
as for the non- adjusted SURF, for a number of different 
z a to generate the adjusted HF field s&(r,t) = s+(r, t) — 

i)z a {s-(r,t-Ta)}. 

4) Fundamental imaging beam: Standard fundamental 
imaging transmit fields without LF manipulation or SURF 
processing. The transmit pulse is equal to the HF part of 
the SURF pulse complex described above. 

III. Results and Discussion 

The propagation effects caused by the modeled body- 
wall, including focus degradation, are further illustrated 
in Fig[2j which shows beam profiles, transmit-beams, and 
axial pulses for (non- adjusted) SURF and fundamental 
beams. Figure [3] shows adjusted SURF beams simulated 
by homogeneous and inhomogeneous propagation with 6 
different choices of suppression depths z a , using both time- 
shift adjustment or filter- adjustment. The corresponding 
beam profiles along the z axis are shown in Fig. SJ 
The time-development of the SURF on-axis pulses as a 
function of depth are displayed in [5] for the same z a set. 
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Fig. 2. Body-wall model impact on the transmit beams illustrated 
through comparison between non- adjusted SURF and fundamen- 
tal imaging transmit fields, for homogeneous and inhomogeneous 
medium. Top pane row: beam profiles in focus. On the left along 
the x axis, on the right along the y axis. Second pane row from 
top: Transmit-beams, from left to right: SH, FH, Sl(zx), Fl(z x), 
Sl(zy), and Fl(zy) with the abbreviations explained in Table III. I 
(zx) and (zy) refer to two planes perpendicular to the beam axis. 
Third pane row from top: Beam-axis pulses, from left to right: SH, 
FH, SI, and FI. Bottom two panes: beam profiles along the depth 
axis. Individual normalization on top and common normalization 

below. Line notations: : fundamental (inhomog.), : 

SURF (inhomog..), : fundamental (homog.), — — — : SURF 

(homog.). 

TABLE II. 

Abbreviations used in the field descriptions of Fig. [2] U 



are about as viable for inhomogeneous as for homogeneous 
media. 

The adjustment abilities in inhomogeneous media are 
further quantified by the specific reverberation suppres- 
sion beam quality ratios Q Za corresponding to the shown 
beams. These are displayed for different suppression 
depths in Fig. [7] using time-shift adjustment and filter- 
adjustment as well as for non- adjusted and fundamental 
beams. This is done both for homogeneous medium and 
for propagation through the body-wall. The figure also 
shows general reverberation suppression quality ratios Q 
for different time-shift adjustments for the same medium 
cases. The specific quality ratio was defined in |2[ as: 



Zf oo 2tt 



OO 2-7T 



(4) 



z=z n r=0 0=0 



=0 0=0 



where E(f) is the beam energy at the spatial coordinate 
r described by the cylindrical coordinates (z,r, 6). The 
imaging depth region is within z G [z n ,Zf] and is equal to 
the focal region with the given aperture and imaging fre- 
quency. The measure Q Za illustrates how well the transmit 
beam is adapted to suppression of reverberations where 
the first scattering takes place at the depth z = z a . The 
general quality ratio was defined in [l| as: 

2/ CX3 2f / Z n CO 2-7T 

«=EEE^/EEE^- (5) 
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These pulses are further inspected in Fig. [6j where their 
signatures at the z a depths are compared to the signatures 
in focus. Side-by-side comparison of the adjusted beams 
for different z a is of particular significance, as it visualizes 
that for the selected z a depths, the adjustment methods 



The measure Q illustrates the ability of the transmit beam 
to suppress reverberations when the first scatterings are 
distributed within the near-field region z G [0, z n ]. A 
related quality measure was applied in [25] for assessment 
of near-field echo suppression in tissue-harmonic imaging. 

The body-wall model chosen causes severe aberration 
disturbance of the transmit fields. Its characteristics, as 
displayed in Table |Lj may be compared to values found 
at body temperature in experimental measurements on 
human abdomen specimens done by Hinkelman et al. (26[, 
where the amplitude correlation length is within 1.3 to 
2.9 mm, the amplitude RMS value within 2.9 to 3.5 dB, 
the time-delay correlation length within 24 to 64 ns, and 
the time-delay correlation length within 3.3 to 17 mm. 
Especially the time-delay RMS fluctuation is more severe 
for the aberrative model utilized here than for the cited 
measurements. 

The RMS magnitude of the time-delay fluctuations 
introduced by the utilized aberrating layer model are, as 
indicated in Table |TJ 450 ns, while the accumulated time- 
delay between s + (r,t) and s_(r,t) needed to generate 
maximum SURF imaging gain following the model intro- 
duced in [l|, is smaller: -^j^ = 142 ns. The SURF synthetic 
transmit-beam is however shown to be generated despite 
these great delays induced by the body-wall. Two factors 
explain this: 

1) The s + and S- propagate following the same path 
and therefore experience the same aberration time-delay 
at each spatial point. Therefore this abberation delay 
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Fig. 3. Adjusted SURF beam cuts for z a = 1.5 cm (top left pane group), z a = 4.5 cm (top right), z a = 10.5 c m (m iddle left), z a = 19.5 cm 
(middle right), z a = 30.5 cm (bottom left), and z a = 39 cm (bottom right). Abbreviation codes shown in Table [TlTl Beams, from left to right 
within each z a pane group: rjSH, rSH, fySI(zx), tSI(zx), fySI(zy), and rSI(zy). The bottom-right colorbar indicates the 35 dB dynamic range 
that is used for all the above beam cuts. 
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Fig. 4. On-axis SURF beam profiles as function of depth z for z a = 1.5 cm (top left pane group), z a = 4.5 cm (top right), z a = 10.5 cm 
(middle left), z a = 19.5 cm (middle right), z a = 30.5 cm (bottom left), and z a = 39 cm (bottom right). Line descriptions: — — — : time-shift 
adjustment (homog.), : time-shift (inhomog.), : filter \) Za adjustment (homog.), : filter \) Za adjustment (inhomog.) 



is canceled when forming the difference field s&{f,i) = 
s+(r,t) - s_(r,t). 

2) The LF part of each SURF pulse complex experience 
the same aberration delay during propagation. Therefore 
the HF pulse remains within roughly the same LF pressure 
as when propagating without aberration. This is illus- 
trated in Fig. [8j where the total propagating SURF pulse 
complexes with positive and negative LF polarities are 
sampled at one depth within the body-wall. 

The SURF transmit-beam in the homogeneous medium 
has suppressed focal sidelobes compared to the fundamen- 
tal transmit-beam, while the mainlobe stays unchanged. 
The suppression is only around 1 dB for the first sidelobe, 
while it increases to around 3 dB for the second sidelobe. 



In the inhomogeneous medium case, the sidelobes are also 
suppressed for SURF compared to fundamental, however 
about half as much as for the homogeneous medium. 

Figure [7] indicates that, both for homogeneous and 
inhomogeneous medium, the time-shift adjustment gives 
superior, more robust and more predictable specific rever- 
beration suppression quality ratio Q Za than the filter-shift 
adjustment for depths where z a >8mm, while the ratio is 
higher for filter- adjustment for shallower depths. The time- 
shift adjusted beams have Q Za around 17 dB above Q Za for 
the non- adjusted SURF beam, both for homogeneous and 
inhomogeneous medium. 

The Q Za ratio calculated in the homogeneous medium 
for the time-shift adjusted beams, the non- adjusted SURF 
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Fig. 5. On-axis SURF pulses as function of depth z and time t for z a = 1.5 cm (top left pane group), z a = 4.5 cm (top middle), z a = 
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Fig. 6. On-axis pressure in focus (bottom sub-panes of each pair) and at z a (top sub-panes of each pair) as function of time t for z a = 1.5 cm 
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Fig. 7. Reverberation suppression beam quality measures with the 
specific ratios Q Za as a function of z a in the top pane and the general 
Q as a function of the adjustment time-shift r a in the bottom pane. 

Time-shift adjusted: (homogeneous) and (inhomoge- 

neous. Filter-adjusted: (homog.) and (inhomog.), No 

adjustment: «■■■■■ (homog.) and i«™ (inhomog.), Fundamental 
beam: — — — (homog.) and (inhomog.). 



beams and the fundamental beams align well with what 
was found for the same 2-D axisymmetric simulations in 
Q. The filter- adjusted SURF beams do however show 
worse and more unstable Q Za in the present 3-D sim- 
ulations than in the referred work. This is due to the 
averaging over the region S which is used here on both 
s+(r, t) and s_(r, t) before \) Za is calculated. Then the 
canceling of temporal parts of the pulses which vary the 
most over the region that is the edge wave parts of 
the pulse, are less suppressed. The amplitude of these 
edge waves may be decreased if apodization is applied to 
the surface of pulse transmission [27], [H|- That way the 
lateral region where \) Za generates high suppression might 
get larger. 

For z a > 8 mm, Q Za is lower both for the homogeneous 
and inhomogeneous media for time-shift, compared to for 
filter-adjustment. Therefore use of time-shift is preferred 
when z a is far from the transducer. This is also reflected 
in the z a = 39 cm section of Fig. [3] both with and without 
the body-wall. There significant in-focus field degradation 
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x [mm] 

Fig. 8. Propagating SURF pulse complexes sampled in the (zx) 
plane at the depth z = 6 mm. Horizontal axis: time, vertical axis: 
lateral direction. Top and middle panes: total SURF pulse complexes 
with opposite LF polarities. Bottom pane: extracted HF difference 
pulse sa(x, t). 

is observed for the filter- adjustment beam as opposed to 
for the time-shift adjusted beam. 

Ultrasound images of patients whose anatomy cause 
noisy images, for example due to obesity, are likely to be 
blurred owing to both aberration and reverberation noise. 

The repetitive nature of reverberation noise makes the 
receive signals correlated both in the temporal and the 
spatial directions. This may obstruct estimation of the 
filter that is used in aberration correction schemes like 
[29( | . By suppression of reverberation noise, estimation 
of the aberration correction filter is thus enhanced, and 
therefore the aberration correction itself is likely to get 
more accurate. Instruments that may combine aberration 
correction with reverberation suppression are thus, com- 
pared to instruments with aberration correction alone, 
likely to produce images of enhanced quality also regarding 
the wave- front aberration noise. 

IV. Concluding remarks 

This work indicates through computer simulations that 
generation a SURF reverberation suppression synthetic 
transmit-beam is feasible in spite of the propagating waves 
being distorted by an inhomogeneous medium which emu- 
lates a strongly aberrating body-wall. Generation of the 
beam is hence also likely to be possible when utilizing 
body-wall models producing less severe aberrations cor- 
responding to what was measured e.g. in [26|. 

The post-processing suppression depth adjustment 
methods, as previously implemented in [2| for a homo- 
geneous medium, are here shown to be feasible on the 
simulated aberrated fields also within an inhomogeneous 
medium. 

Further research of interest within the field also includes 
in vivo imaging through true body-walls using the inves- 
tigated methods and also performance comparisons with 
tissue harmonic imaging, especially when utilizing pulse 
inversion. 

Non-adjusted and adjusted beams could also be gen- 
erated from SURF transmit fields measured by a hy- 
drophone in a water-tank. Beam generation with a in- 
homogeneous body-wall setup may be emulated using 



some aberrating material within the tank, e.g as done 
by use of gel layers in [30] . It is furthermore desirable to 
perform numerical propagation simulations using a tool 
that properly handles multiple scattering (as does not 
the forward-propagation method used here) in order to 
further quantify the reverberation suppression ability of 
the methods within an aberrative medium, as well as 
doing comparisons to the performance of tissue harmonic 
imaging using the same setup. 

We expect the family of SURF reverberation suppres- 
sion methods to come out as versatile techniques to 
enhance ultrasound image reconstruction also in clinical 
settings. 
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